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生 度 对 高 超声 速 钝 锥 边界 层 稳定 性 的 影响 


秦 丽 华 , 韩 宇 峰 , 曾 伟 
(天 津 大 学 力学 系 ,300072 天 津 ) 


摘 要 :采用 直接 数值 模拟 方法 计算 了 8 个 不 同 球 头 半径 的 钝 锥 基本 流 , 运 用 线性 稳定 性 理论 分 析 

了 印度 对 边界 层 稳 定性 的 影响 。 结 果 表 明 , 随 印度 增 大 ,边界 层 内 的 不 稳定 区 向 下 游 移动 ,第 二 模 

人 在 线性 稳定 性 分 析 的 基础 上 ,研究 了 非 线性 扰动 演化 以 及 平均 流 修 正 对 稳 
定性 的 影响 。 结 果 表 明 ,在 基本 流 中 引入 有 限 幅 值 扰 动 后 ,下 游 的 平均 流 剖 面 会 发 生 明 显 改 变 。 流 

场 稳定 性 发 生 显著 变化 ,线性 阶段 最 不 稳定 的 第 二 模 态 波 变 得 稳定 ,而 第 一 模 态 波 明显 增长 起 来 。 

第 一 模 态 波 的 快速 增长 使 V 值 可 以 达到 4, 这 将 会 对 转换 有 很 大 的 促进 作用 。 
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Effect of bluntness on the stability of hypersonic 
blunt cone boundary layers 


QIN Lihua, HAN Yufeng,CAO Wei 


(Department of Mechanics ,Tianjin University ,300072 Tianjin, China) 


Abstract :Instabilities of boundary layers over blunt cones with eight nose radiuses are numerically stud- 
sied. The effect of bluntness on the boundary layer stability is analyzed by linear stability theory. lt is found 
that the passivation of the cone nose leads to a downstream movement of the unstable region and a decrease 
of the maximum growth rate of the second mode. On this basis, nonlinear disturbance evolution and the 
effects of the mean flow distortion on the stability are studied. The results reveal that when the disturbance 
with finite-amplitude is introduced into the flow ,the downstream mean flow profiles are obviously modi- 
fied. As a result, the second-mode, which is the dominant instability in linear growth stage, is damped 
greatly while the first mode wave is amplified significantly. The rapid growth of the first mode can reach 
N=4,which will considerably promote the transition. 

Key words :hypersonic boundary layer; bluntness; direct numerical simulation ;stability analysis; transition 


prediction ; nonlinear disturbance evolution 


边界 层 转换 是 指 边界 层 流动 由 层 流 状态 转变 为 ”” 理 现 象 。 转 所 过程 中 飞行 器 表面 的 摩擦 系数 急剧 升 
济 流 状态 的 过 程 ,是 由 多 个 因素 影响 的 复杂 流动 物 高 ,显著 改变 气动 力 和 壁面 热流 ,直接 影响 飞行 噩 的 
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气动 性 能 和 热 防护 设计 。 因 此 在 高 超声 速 飞行 器 设 
计 中 ,准确 预测 转换 位 置 就 变 得 非常 重要 。 

近年 来 ,从 实验 、 理 论 和 数值 方面 对 转换 进行 了 
大 量 研究 。 如 壁面 温度 对 平板 边界 层 转换 的 影响 以 
及 转换 测量 技术 的 研究 “> ,但 对 于 高 超声 速 钝 锥 边 
界 层 转换 的 研究 相对 较 少 。 如 图 1 所 示 , 给 出 了 钝 
锥 的 示意 图 ,不 同 的 R* 表示 不 同 的 印度 ,R' 越 大 则 
印度 越 大 。 钝 度 是 影响 转 按 的 重要 因素 之 一 ,其 中 ， 
STETSON 等 ?在 马赫 数 5.5 的 实验 中 ,发 现在 小 印 
度 范围 内 增加 钝 度 会 导致 转换 推迟 , 当 钝 度 大 于 某 
个 临界 值 时 ,转换 位 置 反而 提前 ,此 现象 被 称 为 “ 转 
换 反 转 ” 。 后 来 ,SOFTLEY“” 在 实验 中 也 发 现 了 反 
转 现象 。 对 于 小 钝 度 范围 内 转换 推迟 的 现象 可 以 通 
过 由 于 头 部 印度 产生 的 炳 层 使 当地 雷诺 数 减 少 来 解 
释 3 但 对 于 大 钝 度 时 产生 转 按 反 转 的 原因 是 未 知 的 。 
1983 年 ,STETSON 等 ”在 马赫 数 7.99 的 钝 锥 流动 
稳 宕 性 实验 中 未 观察 到 转换 反 转 现象 ,原因 是 实验 
路 合用 的 钝 锥 模型 不 够 长 ,不 足以 发 生 转 按 。 他 们 
的 实验 表明 边界 层 的 扰动 主要 是 第 二 模 态 波 。 最 
他 NATO STO 组 AVT-240 研究 了 钝 度 对 高 超 音速 
角 锥 体 转换 的 影响 ,通过 瞬 态 增长 理论 分 析 发 
珑 向 锥 头 部 附近 的 干扰 会 产生 相对 显著 的 非 模 态 
放大 ,转换 的 开始 与 钝 度 和 粗糙 度 的 综合 影响 
有 se 。 


图 1 钝 锥 示意 氏 
Fig.1 Schematic diagram of blunt cone 

在 数值 模拟 方面 ,关于 钝 度 对 转换 的 影响 研究 
大 都 是 基于 线性 稳定 性 理论 ,依据 STETSON 等 5 马 
赫 数 7. 99 的 实验 工 况 ,许多 学 者 进行 了 数值 研 
究 " 1 。 结 果 表明 , 头 部 钝 化 使 边界 层 更 加 稳定 ,但 
并 未 发 现 转换 反 转 。ZHONG 等 (对 STETSON 
等 马赫 数 5.5 流动 转换 实验 观察 到 的 反 转 现象 进 
行 线性 稳定 性 分 析 。 数 值 研究 了 0. 156、0.5、1.5 英 
才 的 3 个 不 同 头 部 半径 ,但 没有 发 现 反 转 现象 。 由 
于 实验 环境 中 的 噪声 干扰 是 不 可 控 的 ,其 幅 值 可 能 
大 到 使 扰动 增长 绕 过 线性 区 而 迫使 转换 提前 发 生 。 
他 们 提出 有 必要 进一步 研究 非 线性 效应 对 转换 的 
影响 。 

总 之 ,对 于 STETSON 等 写实 验 中 出 现 的 转换 反 
转 现象 ,很 多 学 者 都 进行 了 研究 ,但 都 未 找到 转换 反 
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转发 生 的 原因 。 对 于 钝 锥 , 目前 的 数值 研究 多 集中 
在 线性 稳定 性 分 析 , 而 非 线性 扰动 演化 的 研究 较 少 ， 
因此 ,本 研究 主要 采用 考虑 非 线性 的 扰动 直接 数值 
模拟 方法 ,对 不 同 钝 度 球 锥 的 扰动 演化 进行 分 析 , 以 
研究 非 线性 作用 对 钝 锥 边界 层 稳定 性 的 影响 。 


1 数值 方法 


1.1 物理 模型 


本 研究 的 计算 模型 如 图 2 所 示 ,无 穷 远 来 流 为 
零 攻 角 ,图 中 (x,r) 是 柱 坐 标 系 。x、r 分 别 代表 柱 坐 
标 系 的 轴 癌 与 径 向 ,坐标 原点 选 在 钝 锥 的 球 头 驻 点 
处 。( ,mn ) 是 贴 体 坐 标 系 ,分 别 表 示 沿 锥 面 方向 和 
垂直 锥 面 方向 。0 表示 钝 锥 的 半 锥 角 ,R" 表示 钝 锥 
体 头 部 半径 。 虚 线 包含 的 区 域 为 计算 域 ,点 划 线 表 
示 马 形 激 波 的 大 致 位 置 ,计算 域 的 内 边界 是 钝 锥 体 
的 壁面 ,计算 域 的 外 边界 保证 包含 马 形 激 波 。 
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图 2 物理 模型 的 坐标 系 及 计算 域 示 意图 
Fig.2 The physical model and calculation domain 

图 3 给 出 了 计算 钝 锥 基本 流 的 网 格 分 布 示 例 ， 
流向 和 法 向 网 格 分 布 按 指 数 形 式 分 别 在 圆锥 头 部 和 
壁面 处 加 密 , 网 格 类 型 选择 的 是 光滑 分 布 。 对 于 不 
同 钝 度 的 网 格 分 布 ,根据 其 边界 层 的 厚度 调整 网 格 
的 密度 ,保证 边界 层 内 和 头 部 区 域 有 足够 的 网 格 点 。 

5 Sy 


图 3 计算 网 格 分 布 图 


Fig.3 Distribution of computational grid 
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1.2 计算 方法 及 参数 选取 


本 研究 基本 流 的 计算 采用 柱 坐标 系 下 的 N-S 方 
程 ,由 于 本 研究 的 计算 模型 为 零 攻 角 钝 锥 模型 ,其 流 
动 沿 周 向 一 致 ,在 柱 坐标 系 下 方程 可 以 简化 为 二 维 
形式 ,具体 如 下 。 

OY (1) 

ot Ox or Ox or 
其 中 ;xr 为 轴 向 和 径 向 坐标 ;i 为 时 间 变量 ;U 为 守 
恒 变 量 ;BF 分 别 为 无 条 通 量 ;, ,为 条 性 通 量 ;M 
为 柱 坐 标 系 产生 的 源 项 。 各 项 具体 表达 形式 参考 文 
献 [11] 。 计 算 时 对 方程 式 (1) 中 的 量 进行 无 量 纲 
化 ,其 中 :*v 以 钝 锥 头 部 半径 进行 无 量 纲 化 ;pu 
7 分 别 以 自由 来 流 的 密度 p= 、 速 度 U; 和 温度 7 进 
行 无 量 纲 化 ;压力 p 用 p2 U2 U; 进行 无 量 纲 化 ; 黏 性 
系数 以 迪 进行 无 量 纲 化 ;热传导 系数 由 黏 性 系数 和 
痢 庆 特 数 给 出 ;时 间 ! 和 频率 w 用 有 R /U2 进行 无 量 
纲 化 ;上 标 ” * "表示 有 量 纲 量 ;下 标 * "表示 无 穷 远 
义演 参数 。 

O 〇 数值 计算 时 , 需 将 柱 坐 标 系 下 的 控制 方程 变换 
丽 算 坐标 系 下 。 贴 体 计算 坐标 为 (8,m) ,变换 后 
的 控制 方程 为 

CN 


aF 3 五 aF, 贡 
DE +M 
df dE nm KE 

RU = UE = IEB+EF) ;P= JIE+ 


(2) 


太 本 ) ;M = JM; J 为 坐标 变换 的 雅克 比 行列 式 。 
“基本 流 的 计算 中 ,壁面 边界 条 件 采 用 无 滑 移 和 
等 温 壁面 条 件 , 极 轴 处 使 用 轴 对 称 边界 条 件 ,上 边界 
给 定 自 由 来 流 条 件 ,出 口 采用 基本 量 线性 外 推 条 件 。 
扰动 演化 的 计算 是 在 基本 流 的 基础 上 引入 有 限 幅 值 
的 扰动 ,对 于 二 维 情况 ,入 口 边 界 条 件 为 
u(t,0,y) = uo(0,y) + A [a(y)e 7] 
v(t,0,y) 过 20(0,yY) +4 [2(y)e 2 ] 


T(1,0,7) = T,(0,y) + ATT(y) ee 0],, 


p(t,0,y) = po(0,y) + ATP(Yy)e ， 


(3) 
其 中 : L(y) 、 人 (y) 、7T(y) 、P(y) 分 别 为 对 应 物理 
量 的 特征 函数 ; w 为 人 口 给 定 的 频率 ;4 为 扰动 的 幅 
值 ;下 标 real 表示 取 实 部 。 出 口 处 使 用 骨 边 区 条 
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件 ' 中 ,将 扰动 在 出 口 处 吸收 衰减 为 零 ,避免 扰动 反 
射 对 上 游 流 场 产 生 影响 。 由 于 本 研究 计算 的 扰动 是 
二 维 平面 波 , 展 向 波 数 B =0, 扰 动 波 在 向 下 游 演 化 
的 过 程 中 始终 是 沿 圆锥 母线 传播 的 平面 波 。 柱 坐标 
系 下 的 N-S 方程 可 以 简化 为 二 维 N-S 方程 ,因此 , 扰 
动 演 化 的 计算 也 采用 二 维 N-S 方程 。 

计算 参数 选取 与 STETSON 等 中 实验 工 况 相同 ， 
零 攻 角 , 半 锥 角 8 度 , 来 流 马 赫 数 Ma, =5.5 ,Pr = 
0.72,y = 1.4, 壁面 温度 7, = 296 K, 其 他 参数 见 
表 1。 计 算 时 对 流 项 应 用 五 阶 WENO 格式 , 黏 性 项 使 
用 六 阶 中 心 差分 格式 , 时间 项 采用 三 阶 Runge- 
Kutta 法 。 


表 1 自由 来 流 参数 表 


Tab.1 Computational parameters 


密度 / 温度 / 压力 / 动力 黏度 / ”单位 雷诺 数 / 
(kg .m-3) K Pa (N.s.m-2) (1.m-!) 
1.549 x107! 174.46 7756.56 1.176x10- 1.895 x107 


1.3 ” 转 拨 预 测 方法 


对 于 转换 位 置 预测 ,本 研究 采用 的 是 基于 线性 
稳定 性 理论 的 e 方 法” ,根据 线性 稳定 性 理论 ， 
将 小 扰动 写成 行进 波 的 形式 

P(E MLst) = P(N) ee (3) 

其 中 :上 7 分别 代表 流向 坐标 法 向 坐标 、 周 向 坐 

标 ;w 是 扰动 波 的 频率 ; a 和 BB 分别 为 流向 波 数 和 展 

向 波 数 ,对 于 空间 模式 a 为 复数 , -a 表示 扰动 的 增 

长 率 。 对 增长 率 沿 流向 积分 可 得 到 幅 值 放大 倍数 
N,N 值 表达 式 为 

N= (元 ) = 一 | ad (4) 


其 中 :46 为 扰动 的 初始 幅 值 ;4 为 扰动 增长 的 幅 值 。 
对 于 高 超声 速 边界 层 转 按 预测 现在 还 没有 具体 的 实 
验 数据 能 够 确定 N 值 取 多 大 作为 转换 判 据 。 本 研究 
依据 改进 的 ev 方 法 , 认为 扰动 的 初始 幅 值 都 为 
0.3% , 当 幅 值 增 大 到 1. 5% 时 认为 转换 开始 ,此 时 ， 
根据 式 (4) 可 以 计算 出 N = Im( 祭 )= 4。 


1.4 程序 验证 


基本 流 计算 程序 和 线性 稳定 性 理论 程序 的 可 靠 
性 验证 ,选择 与 文献 110] 相同 的 工 况 ,计算 了 4、 
12.7.38.1 mm,3 个 印度 的 钝 锥 基本 流 场 ,并 进行 稳 
定性 分 析 。 图 4 给 出 了 钝 度 为 4 mm 的 基本 流 的 流 
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向 速度 剖面 对 比 ,图 5 给 出 了 3 种 工 况 的 中 性 曲线 
对 比 图 。 可 以 看 出 与 文献 [10] 的 计算 结果 吻合 较 
好 ,表明 本 研究 的 程序 是 可 靠 的 。 
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流向 速度 剖面 与 文献 [10] 的 结果 对 比 (s” =0.4 m) 
Fig.4 Comparison of flow velocity profile 
with Ref. [10|] (s* =0.4 m) 
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图 5 中 性 曲线 与 文献 [10] 的 结果 对 比 

Fig.5 Comparison of neutral curves with Ref. [10] 

对 于 扰动 计算 程序 的 验证 ,由 于 小 扰动 直接 数 
值 模拟 ( DNS) 的 结果 应 该 与 线性 稳定 性 理论 (LST) 
的 结果 相符 。 为 此 ,选择 印度 4 mm 的 情况 ,从 基本 
流 取 140 倍 到 210 倍 的 头 部 半径 区 域 作 为 计算 域 。 
在 计算 域 入 口 引 入 基本 流 线 性 稳定 性 分 析 得 到 的 第 
二 模 态 波 ,其 频率 为 w = 11 ,初始 幅 值 为 1 x10…， 
计算 到 周期 定常 后 得 到 扰动 幅 值 增长 曲线 ,并 与 线性 
稳定 性 理论 预测 的 幅 值 增长 结果 对 比 。 如 图 6 所 示 ， 
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实 线 表示 DNS 的 结果 ,虚线 表示 LST 预测 的 结果 ,可 
以 看 出 两 者 吻合 较 好 ,表明 扰动 计算 程序 可 靠 。 


“~ 
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3 


图 6 DNS 与 LST 幅 值 曲线 对 比 图 
Fig.6 Comparison of amplitude between DNS and LST 


2 基本 流 计算 结果 与 分 析 


在 先前 的 研究 中 ,大 多 数 数值 计算 只 选择 了 
3 个 钝 度 (LEI 等 " ,2012 ) ,但 都 未 发 现 转换 反 转 。 
因此 ,在 实验 工 况 的 基础 上 ,本 研究 进一步 扩大 了 钝 
度 范围 。 选 取 了 8 个 不 同 钝 度 ,分 别 为 4、12.7、20、 
30 38.1 .50 .60 80 mm, 其 中 4、12.7、38.1 mm 为 实 
验 中 的 3 个 印度 。 对 不 同 印度 的 基本 流 本 研究 都 进 
行 了 网 格 无 关 性 验证 ,法 向 网 格 边 界 层 内 不 少 于 
101 个 网 格 点 ,流向 网 格 头 部 区 域 做 加 密 处 理 , 保 证 
计算 结果 的 可 靠 性 。 由 于 扰动 演化 的 研究 需要 在 线 
性 稳定 性 分 析 的 基础 上 进一步 开展 ,因此 本 章 首 先 
介绍 线性 稳定 性 分 析 的 一 些 计 算 结 果 , 为 进一步 数 
值 模拟 扰动 演化 做 准备 。 下 面 给 出 了 不 同 钝 度 情况 
下 的 基本 流 和 稳定 性 分 析 结 果 的 对 比 。 


2.1 无 量 纲 坐 标 下 基本 流 场 特 性 及 稳定 性 的 分 析 


由 于 无 量 纲 特征 长 度 选取 的 是 圆锥 前 缘 头 部 半 
es 

径 ,根据 公式 Re = “一 可 知 ,在 无 量 纲 坐标 系 

下 , 钝 度 的 变化 通过 改变 来 流 雷 诺 数 来 实现 , 钝 度 对 

基本 流 场 的 影响 体现 为 雷诺 数 变化 对 流 场 特 性 的 影 

响 。 不 同 钝 度 下 的 来 流 雷诺 数 见 表 2,R" 表示 半径 ， 

Re 。 表示 来 流 雷诺 数 。 


表 2 不 同 印度 的 来 流 雷 诺 数 


Tab.2 Free-stream Reynolds number for different bluntness 


站 专 4 mm 


12.7 mm 20 mm 30 mm 


38.1 mm 50 mm 60 mm 80 mm 


Re, 75 800 240 665 379 000 


图 7 给 出 了 沿 流向 同一 位 置 = 150(s 表示 从 


568 500 


721 995 


圆锥 头 部 驻 点 起 沿 锥 面 方向 的 坐标 ,下 同 ) ,不 同 钝 


947 500 1 137 000 1516 000 
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度 的 钝 锥 基本 流 的 流向 速度 和 温度 剖面 。 由 图 
7(a) 可 以 看 出 随 着 钝 度 的 增加 ,速度 剖面 越 来 越 靠 
近 壁 面 ,表明 边界 层 变 得 越 来 越 薄 。 图 7(p) 也 可 看 
出 , 随 着 印度 的 增加 边界 层 变 薄 , 同 时 边界 层 内 的 温 
度 剖 面 峰值 变 大 。 由 于 头 部 半径 的 增 大 ,雷诺 数 也 
随 之 增 大 ,流体 黏 性 力 变 小 使 得 边界 层 内 的 黏 性 区 
越 来 越 小 ,边界 层 变 薄 。 


0.3 


一 -一 4mm 


一 一 12.7 mm 


(a) 速度 剖面 


一 一 4mm 


一 2.7imm 
一 一 一 一 一 20 mm 
一 一 30mm 
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(b) 温度 剖面 


图 7 无 量 纲 坐标 下 不 同 钝 度 基 本 流 的 速度 剖面 
和 温度 剖面 (* =150) 


Fig.7 Velocity and temperature profiles in nondimensional 


coordinates with different bluntness (s =150) 

根据 线性 稳定 性 理论 ,对 不 同 钝 度 的 定常 基本 
流 进行 稳定 性 分 析 , 计 算得 到 不 同 钝 度 下 基本 流 中 
性 曲线 和 最 大 增长 率 治 流向 的 变化 ,如 图 8 所 示 。 

图 8(a) 为 不 同 钝 度 下 的 中 性 曲线 ,可 以 看 出 不 
同 钝 度 的 基本 流出 现 不 稳定 波 的 频率 区 间 有 很 大 的 
差别 , 钝 度 最 小 的 不 稳定 波 频 率 的 范围 也 最 小 , 随 着 
印度 增 大 ,不 稳定 波 出 现 的 范围 随 之 变 大 。 从 整体 
上 看 ,印度 越 大 ,最 先 出 现 不 稳定 波 的 位 置 也 跟着 问 
靠近 钝 锥 头 部 的 方向 移动 ,但 移动 的 量 不 是 很 大 。 
图 8(b) 为 最 大 增长 率 沿 流向 的 变化 ,可 以 看 出 , 随 
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钝 度 增 大 ,最 大 增长 率 也 随 之 增 大 。 在 无 量 纲 坐 标 
下 ,不 同 钝 度 的 钝 锥 之 间 存 在 相似 性 。 


120 4 mm 
12.7 mm 
一 一 一 20mm 
100 一 一 .一 .一 一. 30 mm 


(a) 中 性 曲线 


0.7 一 一 4mm 
一 一 一 12.7mm 


(b) 流向 最 大 增长 率 
图 8 不 同 钝 度 基本 流 的 中 性 曲线 和 
最 大 增长 率 沿 流向 的 变化 


Fig.8 Neutral curves and maximum growth rate along the 


streamwise direction for basic flow with different bluntness 
2.2 有 量 纲 坐 标 下 基本 流 场 特性 及 稳定 性 的 分 析 


图 9 给 出 了 有 量 纲 坐 标 下 同一 位 置 不 同 钝 度 基 
本 流 剖 面 的 比较 结果 。 可 以 看 到 随 着 钝 度 的 增加 ， 
有 量 纲 的 边界 层 厚 度 变 大 ,这 是 由 于 钝 度 增 大 使 流 
体 动能 亏损 变 多 ,边界 层 内 速度 变化 到 来 流速 度 的 
趋势 变 得 越 来 越 缓 。 从 温度 剖面 可 以 明显 看 出 , 随 
着 钝 度 增 大 , 激 波 出 现 的 位 置 逐渐 远离 壁面 。 

图 10 给 出 的 是 有 量 纲 坐标 下 不 同 钝 度 基 本 流 
的 中 性 曲线 和 最 大 增长 率 沿 流向 的 变化 。 从 图 
10(a) 可 以 更 直观 地 看 到 钝 度 对 边界 层 稳定 性 的 影 
响 , 钝 度 越 大 ,基本 流 的 不 稳定 区 域 越 向 下 游 移 动 ， 
不 稳定 波 的 有 量 纲 频率 最 大 值 也 越 来 越 小 。 图 
10(b) 给 出 了 有 量 纲 最 大 增长 率 随 流向 的 变化 , 随 
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着 印度 增 加 ,有 量 纲 最 大 增长 率 明 显 减 小 。 这 都 表 
明 第 二 模 态 波 随 钝 度 增 大 越 来 越 稳定 , 钝 度 对 第 二 


模 态 起 稳定 性 作用 。 
10 
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(b) 不 同 钝 度 的 流向 温度 剖面 

合 9 有 量 纲 坐 标 下 的 速度 剖面 和 温度 前 面 (\” =2 m) 

Shig. 9 Velocity and temperature profiles in dimensional 
© coordinates with different bluntness (s* =2 m) 

站 图 11 给 出 了 不 同 印度 基 本 流通 过 线性 稳定 性 
分 析 得 到 的 图 值 曲线 。 可 以 看 出 ,NN 值 曲 线 的 斜率 
随 钝 度 的 增 大 而 减 小 ,表明 取 同 一 w 值 作为 转 拨 判 
据 时 ,转换 位 置 随 钝 度 增加 而 向 下 游 移动 。 
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(b) 最 大 增长 率 
图 10 有 量 纲 坐 标 下 中 性 曲线 和 最 大 增长 率 沿 流向 的 变化 


Fig. 10 ”Neutral curves and maximum growth rate along 


the streamwise direction in dimensional coordinates 
30 
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人 2.7 Mm 
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50 mm 


sm 
图 11 不 同 钝 度 基本 流 的 V 值 曲线 
Fig.11 N-factor curves of basic flow with different bluntness 
在 图 12 中 给 出 了 不 同 钝 度 下 基本 流 的 转换 位 
置 与 头 部 半径 的 关系 ,可 以 发 现 V 值 分 别 取 4、5、6 
和 7 时 得 到 的 规律 是 相似 的 。 即 转换 位 置 随 钝 度 的 
增 大 而 向 下 游 移动 。 
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图 12 ”转换 位 置 与 头 部 半径 关系 曲线 


Fig. 12 The transition position at different N-factors 
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在 文献 [10] 计算 的 3 个 钝 度 工 况 中 (参照 图 
5) ,得 到 的 稳定 性 分 析 结果 表明 ,第 二 模 态 不 稳定 波 
在 基本 流 中 最 先 出 现 的 位 置 随 印度 的 增 大 向 下 游 移 
动 。 在 本 研究 的 8 个 不 同 钝 度 的 钝 锥 基本 流 中 , 尽 
管 扩 大 了 钝 度 范 围 ,但 通过 线性 稳定 性 分 析 得 到 的 
结果 与 现 有 的 数值 研究 结果 一 致 ”” ,圆锥 前 缘 钝 化 
使 得 转换 推迟 ,在 数值 计算 中 尚未 获得 "转换 反 转 ” 
现象 。 为 了 进一步 研究 印度 对 边界 层 稳 定性 的 影 
啊 , 在 下 一 节 将 介绍 非 线性 扰动 对 稳定 性 的 影响 。 


3 ”有 限 幅 值 扰动 演化 的 数值 研究 


对 于 小 扰动 波 的 扰动 演化 ,其 对 基本 流 场 不 会 
产生 影响 , 流 场 的 流动 稳定 性 与 基本 流 相同 ,但 由 于 
实际 实验 环境 的 复杂 性 ,扰动 可 能 存在 着 多 样 性 和 


流 的 基础 上 进一步 开展 了 有 限 幅 值 扰动 演化 的 数值 
300 a 
拟态 动 演化 计算 结果 


局 计 算 域 的 选取 是 在 原来 的 钝 锥 基本 流 计算 域 基 
础 汪 重 新 截取 一 段 计 算 域 ,作为 扰动 演化 的 新 计算 
达 新 的 计算 域 如 图 13 所 示 , 以 此 计算 域内 的 流 场 
汶 基 本 流 ,计算 域 和 人口 距离 钝 头 前 缘 约 140 个 头 
部 盏 径 。 
EF EN 人 
[® pi By TO LULZ 


New-domain 
Original domain 
Shock 
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图 13 有 限 幅 值 扰动 演化 的 计算 域 
Fig. 13 The computational domain for 
the finite-amplitude disturbance 
以 头 部 半径 4 mm 的 工 况 展示 有 限 幅 值 扰动 演 
化 的 计算 结果 。 在 计算 域 人 口 处 引入 线性 稳定 性 分 
析 得 到 的 最 不 稳定 的 第 二 模 态 波 。 图 14 是 由 线性 
稳定 性 理论 得 到 的 第 二 模 态 特征 函数 曲线 图 ,特征 
函数 以 边界 层 内 速度 特征 函数 幅 值 的 最 大 值 进行 归 
一 化 ,其 实 部 和 虚 部 如 图 14(a) 和 14(b) 所 示 。 


ChinaXiv 合 作 期 刊 第 40 郑 


0.6 0.6r 
0.5 0.5 多 
从 
0.4 人 | 多 
0.3 < 0.3 
0.2 (02 NE 
0.1 1 | 
0.0L -A 0.0 - 
-02 0 02 04 0.6 0.8 1.0 1 0 
AN AN 八 八 
入， 从 Dy 
(a) 速度 特征 函数 (b) 温度 特征 函数 


图 14 入 口 加 入 的 第 二 模 态 特征 函数 

Fig. 14 The second mode eigenfunction introduced at the inlet 

首先 ,研究 了 频率 为 w =11 初始 幅 值 为 1x10” 
的 扰动 。 图 15 给 出 了 扰动 演化 的 结果 ,其 中 图 15 
(a) 是 LST 预测 的 扰动 幅 值 增长 曲线 与 DNS 幅 值 增 
长 结果 的 对 比 ,可 以 看 到 在 x = 162 之 前 DNS 的 结 
果 可 与 LST 的 结果 大 致 重合 ,但 在 x = 162 之 后 DNS 
的 计算 结果 开始 偏离 LST 预测 的 幅 值 增长 趋势 , 且 
DNS 的 结果 明显 小 于 LST 的 结果 。 表 明 此 时 边界 层 
内 已 经 开始 出 现 不 可 忽略 的 非 线 性 扰动 。DNS 的 扰 
动 幅 值 由 1 x10” 增 长 到 0.3 左右 后 开始 趋 于 饱和 。 
从 扰动 的 瞬时 速度 增长 图 15(b) 中 可 以 更 直观 地 看 
到 ,在 x =162 之 后 瞬时 速度 的 增长 开始 出 现 上 下 不 
对 称 ,这 表明 基本 流 开始 出 现 修 正 , 在 下 游 x = 170 
之 后 瞬时 速度 增长 逐渐 趋 于 饱和 。 
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(a) DNS ST 的 扰动 由 (b) 扰动 瞬时 速度 脉动 
图 15 初始 幅 值 为 1 x10 习 的 扰动 演化 计算 结果 
Fig.15 Calculation results of disturbance evolution 

with the initial amplitude of 1 x10 
图 16 给 出 了 初始 幅 值 为 1 x 10 一 的 扰动 演化 结 

果 , 从 中 可 以 发 现 DNS 的 幅 值 增长 结果 更 早 的 偏离 

LST 预测 的 幅 值 增长 曲线 ,大 约 在 x = 150 处 就 发 生 

偏离 ,上 且 在 大 约 x = 160 左右 DNS 的 幅 值 增长 就 趋 

于 饱和 。 还 可 以 看 出 扰动 瞬时 速度 的 增长 相 比 入 口 

幅 值 1 x10 一 的 增长 更 快 ,这 都 表明 在 初始 幅 值 更 大 

的 情况 下 , 非 线 性 扰动 出 现 的 位 置 更 靠 上 游 。 因 此 ， 

为 了 分 析 非 线性 扰动 对 稳定 性 的 影响 ,在 计算 域 人 人 

口 选择 引入 的 扰动 初始 幅 值 为 1 x10。 
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(a) DNS 与 LST 的 扰动 幅 值 增长 对 比 (b) 扰动 瞬时 速度 脉动 
16 初始 幅 值 为 1 x10 习 的 扰动 演化 计算 结 


Fig.16 Calculation results of disturbance evolution 
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with the initial amplitude of 1 x10 
3.2 平均 流 计算 结果 与 分 析 


平均 流 修正 与 转 扎 密切 相关 ,基于 以 上 认 
识 ,首先 对 平均 流 场 进 行 分 析 。 将 计算 到 周期 定常 
的 瞬时 扰动 流 场 对 时 间 平 均 , 以 获得 平均 流 场 。 以 

的 工 况 为 例 , 入 口 扰动 频率 为 w =11, 初 始 扰 

车 值 为 1 x10 习 ,图 17 给 出 了 沿 流向 不 同 站 位 的 
可 以 看 到 , 随 着 扰动 波 向 下 游 

断 演 化 ,流向 速度 剖面 和 温度 剖面 都 出 现 了 很 

国 。 
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(b) 平均 流 温 度 剖面 


图 17 不 同 流向 站 位 平均 流 剖 国 


Fig.17 Mean flow profiles at different streamwise locations 


投稿 网 站 :http://cjam. xjtu. edu. cn 


秦 丽 华 ， 等 : 钝 度 对 高 超声 速 钝 锥 边界 层 稳定 


性 的 影响 hinaXiv 合 作 期 刊 ”27 


图 18 给 出 了 沿 流向 不 同位 置 ,平均 流 剖 面 和 基 
本 流 谢 面 的 广义 拐点 (GIP) 的 分 布 。 可 以 看 出 , 基 
本 流 的 边界 层 中 没有 广义 拐点 ,而 对 于 平均 流 , 在 下 
游 *=180 和 s=200 的 流向 位 置 边界 层 内 出 现 了 4 
个 广义 拐点 。 表 明 在 平均 流 场 的 下 游 其 稳定 性 发 生 
了 明显 改变 。 
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Fig. 18 Generalized inflection points of mean flow 


at different streamwise locations 
运用 线性 稳定 性 理论 对 得 到 的 平均 流 进行 线性 
稳定 性 分 析 , 图 19 给 出 了 沿 流 向 不 同位 置 ,稳定 性 
分 析 后 得 到 的 平均 流 中 性 曲线 (B-w 平面 ) 与 原 基 本 
流 中 性 曲线 的 对 比 。 
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图 19 ”平均 流 与 基本 流 不 同 流向 位 置 的 中 性 曲线 对 比 


Fig.19 Neutral curves of the mean flow and the basic flow 


CC 
人 


at different streamwise locations 
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可 以 看 到 随 着 扰动 向 下 游 的 演化 ,平均 流 的 失 
稳 区 域 逐 渐 偏 离 原来 的 不 稳定 区 。 在 流向 s = 152 
的 位 置 ,平均 流 的 中 性 曲线 在 w = 12 附近 开始 出 现 
凹陷 ,继续 向 下 游 演化 ,在 流向 * = 158 的 位 置 ,平均 
流 的 中 性 曲线 分 为 两 个 不 稳定 区 域 (分 别 用 a、b 表 
示 )。 再 继续 向 下 游 ,除了 高 频 的 两 个 区 域 ,平均 流 
在 低频 位 置 出 现 了 一 个 新 的 不 稳定 区 域 。 

图 20 给 出 的 是 平均 流 不 同 流向 位 置 的 中 性 曲 
线 对 比 图 ,可 以 看 到 平均 流 的 3 个 失 稳 区 域 在 向 下 
游 的 演化 过 程 中 发 生 的 变化 ,其 中 a 和 上 b 两 个 不 稳 
定 区 域 的 范围 都 在 缩小 ,同时 展 向 波 数 也 在 减 小 。 
而 c 区 域 对 应 的 低频 区 的 不 稳定 波 范 围 却 在 扩大 ， 
且 不 稳定 波 的 频率 和 展 向 波 数 都 在 不 断 变 大 。 表 明 
边界 层 内 的 扰动 非 线性 作用 使 得 两 个 高 频 区 的 不 稳 
定 波 变 得 越 来 越 稳定 ,而 低频 区 的 不 稳定 波 增 长 
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i Fig.20 Neutral curves of mean flow 


© at different streamwise locations 

为 了 分 析 这 3 个 失 稳 区 域 扰动 的 特征 ,给 出 了 
s=175 位 置 处 ,平均 流 3 个 失 稳 区 域 对 应 不 稳定 波 
相 速 度 随 频率 的 变化 ,如 图 21 所 示 。 从 相 速 度 可 以 
看 出 频率 最 低 的 不 稳定 波 类 似 于 第 一 模 态 ,而 男 外 
两 个 高 频 区 的 不 稳定 波 类 似 于 第 二 模 态 。 图 22 给 
出 的 是 同一 流向 位 置 ,s =175 处 的 3 个 模 态 对 应 的 
特征 函数 。 可 以 看 出 频率 最 低 的 模 态 对 应 的 速度 和 
温度 特征 函数 峰值 都 远离 壁面 ,这 与 第 一 模 态 特征 
函数 的 特征 相同 。 因 此 ,频率 最 低 的 模 态 对 应 于 第 
一 模 态 。 另 外 两 个 模 态 的 速度 和 温度 特征 函数 最 大 
值 都 靠近 壁面 ,这 与 第 二 模 态 特征 函数 的 特征 相同 ， 
因此 这 两 个 模 态 都 属于 第 二 模 态 ,但 这 两 个 模 态 的 
温度 特征 函数 的 峰值 存在 一 定 的 差别 ,是 两 个 不 完 
全 相同 的 模 态 (分 别 用 第 二 模 态 _a 和 第 二 模 态 _b 表 
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示 )。 表 明 由 于 非 线 性 作用 导致 的 平均 流 修正 ,使 得 
基本 流 中 的 第 二 模 态 分 成 了 两 个 模 态 。 


1.05 


图 21 平均 流 第 一 模 态 和 第 二 模 态 相 速 度 随 频 率 的 变化 
Fig.21 Variation of the phase speed of the first 


and second modes along frequency 


0.0 
0.0 02 04 0.6 0.8 1.0 


向 
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图 22 平均 流 第 一 模 态 和 第 二 模 态 特征 函数 曲线 
Fig.22 Shape-functions of the first and second 


modes in mean flow 

图 23 给 出 了 沿 流向 不 同 站 位 ,3 个 失 稳 区 域 第 
一 模 态 和 第 二 模 态 增长 率 随 频率 的 变化 。 可 以 发 
现 , 随 着 向 下 游 的 演化 第 二 模 态 的 增长 率 不 断 缩 小 ， 
增长 率 最 大 值 从 0. 12 逐渐 缩小 到 0. 06 ,而 第 一 模 态 
的 增长 率 明 显 增 大 ,增长 率 的 最 大 值 从 最 初 的 0.03 
增 大 到 0. 06 左右 。 这 表明 扰动 向 下 游 演化 的 过 程 
中 ,由 于 非 线 性 作用 对 平均 流 剖 面 的 修正 ,使 得 第 二 
模 态 波 的 作用 在 减 小 ,而 第 一 模 态 不 稳定 波 的 作用 
在 不 断 增 大 。 

对 于 其 他 钝 度 , 本 研究 也 进行 了 与 上 述 同 样 的 
计算 ,选取 相同 的 计算 域 人 口 位 置 , 距 前 缘 约 140 倍 
的 头 部 半径 ,在 入 口 处 引入 该 位 置 所 对 应 的 最 不 稳 
定 的 第 二 模 态 波 ,初始 幅 值 相 同 均 为 1 x10 一 。 通 过 
直接 数值 模拟 计算 得 到 平均 流 ,并 对 平均 流 进行 稳 
定性 分 析 。 图 24 给 出 了 12.7、38.1、60 mm,3 个 不 
同 印度 下 的 平均 流 的 第 一 模 态 和 第 二 模 态 增长 率 。 
可 以 看 出 ,对 于 不 同 印度 的 情况 ,出 现 的 影响 规律 类 
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似 , 随 着 扰动 向 下 游 的 演化 ,都 出 现 了 第 一 模 态 增 
大 ,第 二 模 态 减 小 的 现象 。 
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(b) 第 一 模 态 增 长 率 


平均 流 不 同 流向 位 置 第 一 模 态 第 二 模 态 增长 率 曲线 


Fig. 23 The first and second mode growth rate 


na 淮 v 


of the mean flow at different streamwise locations 
斧 从 第 一 模 态 和 第 二 模 态 的 增长 率 曲线 可 以 看 
出 6 第 一 模 态 的 增长 率 已 经 达到 与 第 二 模 态 增长 率 
相同 的 数量 级 。 为 了 分 析 第 一 模 态 对 转 按 的 影响 ， 
图 25 给 出 了 8 个 不 同 钝 度 下 的 第 一 模 态 N 值 曲 线 。 
可 以 看 出 ,对 于 4 mm 的 情况 ,N 值 只 能 达到 2 ,但 对 
于 其 他 钝 度 的 情况 ,N 值 均 大 于 4。 这 表明 有 限 幅 值 
扰动 演化 的 非 线 性 作用 使 得 第 一 模 态 变 得 不 稳定 ， 
可 能 在 转换 过 程 中 起 着 重要 作用 。 

在 基本 流 中 引入 有 限 幅 值 扰动 后 ,基本 流 剖 面 
出 现 明显 修正 , 流 场 不 稳定 性 发 生 了 显著 变化 。 结 
果 表 明 , 边 界 层 中 的 第 二 模 态 变 得 越 来 越 稳定 ,而 第 
一 模 态 的 不 稳定 性 明显 增加 。 对 于 不 同 钝 度 的 情况 
均 存 在 上 述 现象 。 由 于 有 限 幅 值 扰动 的 影响 与 计算 
域 和 人口 位 置 的 选取 以 及 扰动 的 初始 幅 值 均 相关 。 因 
此 ,很 难 通过 e* 方 法 来 定量 的 计算 预测 平均 流 的 转 
按 位 置 。 
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24 不 同 钝 度 下 的 平均 流 的 第 一 模 态 和 第 二 模 态 增长 率 


Fig.24 The first and second mode growth rate of the mean 


flow for different bluntness 
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图 25 不 同 钝 度 下 的 平均 流 第 一 模 态 NN 值 曲线 
Fig.25 The first mode N-factors for different bluntness 


4 结 论 


本 研究 通过 直接 数值 模拟 方法 计算 了 来 流 参数 
与 STETSON 等 "| 实验 工 况 相同 的 8 个 不 同 钝 度 的 
印 锥 基本 流 场 ,分 析 了 印度 对 基本 流 场 特性 的 影响 ， 
并 通过 线性 稳定 性 理论 对 流 场 进 行 了 稳定 性 分 析 。 
在 此 基础 上 研究 了 考虑 非 线 性 扰动 产生 的 基本 流 修 
正 对 流 场 及 稳定 性 的 影响 ,得 到 如 下 结论 。 
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1) 在 4 ~80 mm 的 钝 度 范围 内 ,基本 流 的 线性 
稳定 性 分 析 结 果 表明 , 随 着 钝 度 的 增加 ,边界 层 内 第 
二 模 态 不 稳定 波 的 频率 降低 。 且 边界 层 内 的 不 稳定 
区 域 逐 渐 向 下 游 移动 ,最 大 增长 率 明显 降低 。 通 过 
e 方法 预测 转换 位 置 的 结果 表明 , 随 钝 度 增 大 ,转换 
位 置 向 下 游 移动 ,未 发 现 转换 反 转 现象 。 

2) 引 入 有 限 振 幅 扰 动 后 ,平均 流 剖 面 得 到 明显 
的 修正 。 非 线性 作用 使 得 边界 层 内 的 第 二 横 态 波 越 
来 越 稳定 ,而 第 一 模 态 不 稳定 波 明 显 增长 起 来 。 第 
一 模 态 波 的 快速 增长 使 入 值 可 以 达到 N=4, 这 将 会 
对 转换 有 很 大 的 促进 作用 。 这 与 “转换 反 转 ”的 发 
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